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Abstract. Narrow-Line Seyfert 1 galaxies (NLSls) and Narrow-Line quasars (NLQs) seem to amount to ~ 10-30% of active 
galactic nuclei (AGNs) in the local universe. Together with their average accretion rate, we argue that a black hole (BH) growth 
by factor of 8 - 800 happens in these super-Eddington accretion phase of AGNs. Moreover, there is a possible, systematic 
underestimation of accretion rates (in the Eddington unit) due to an overestimation of BH mass by massive accretion discs for 
super-Eddington objects. If it is true, the factor of BH growth above may be larger by order(s) of magnitude. In contrast, the 
growth factor expected in sub-Eddington phase is only ~ 2. Therefore, the cosmic BH growth by accretion is likely dominated 
by super-Eddington phase, rather than sub-Eddington phase which is the majority among AGNs. 

This analysis is based on the fraction and the average accretion rate of NLSls and NLQs obtained for z ^ 0.5. If those numbers 
are larger at higher redshift (where BHs were probably less grown), super-Eddington accretion would be even more important 
in the context of cosmic BH growth history. 
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1. Introduction 

Almost all galaxies in the local universe have supermassive 
black holes (BHs) in their nuclei, not only active galactic nu- 
clei (AGNs) but also apparently normal galaxies (Kormendy 
& Richstone 1995; Magorrian et al. 1998). Luminous quasars 
with massive BHs are observed till redshift z ~ 6, when the 
universe was less than 10% of its recent age (e.g. Fan et al. 
2003; Vestergaard 2004). However, it is not clear at all when 
those BHs have been formed and how their BH mass (Mbh) 
have evolved. 

Most of unobscured AGNs show broad Balmer lines of hy- 
drogen (FWHM of Hy6 > 2000 km s~'), and they are named as 
broad line Seyfert 1 galaxies (BLSls) and quasars. From sev- 
eral independent arguments, the duration of an AGN is thought 
to be ~ 10 8 yr (Martini 2004 for a review). And the average ac- 
cretion rate M of BLSls/quasars is about 3 L^alc 1 (e.g., Fig. 
1 1 of Kawaguchi 2003), where Lem is the Eddington lumi- 
nosity [~ 10 46 (M BH /10 8 M o ) erg s -1 ]. Thus, BHs increase their 
mass during ~ 10 8 year only by a factor of ~ 2. 

Then, what else can be responsible for the cosmic 
BH growth? Semi-analytical approaches (e.g. KaufFmann & 
Haehnelt 2000) employ a scheme where a BH mass is increased 
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not only by accretion but also by merging of two massive BHs 
during a merger event of galaxies, where the amount of BH 
masses is conserved. However, the BH growth via merging is 
only expected at high mass BHs (associated in most massive 
dark halos; Cattaneo 2002), and thus this mechanism is not 
feasible for BH growth from a seed BH to a Seyfert-level BH 
(~ 10 6 ~ 8 M o ). Moreover, discussions comparing the evolution 
of luminosity functions (LFs) with the volume mass density 
of BHs (Soltan 1982; Chokshi & Turner 1992) infer that the 
BH growth by accretion process is enough to explain the lo- 
cal BH mass density (Fabian & Iwasawa 1999; Yu & Tremaine 
2002; Elvis, Risaliti & Zamorani 2002). Namely, merging of 
BHs might be negligible for the cosmic BH growth (Marconi 
et al. 2004). 

There is growing evidence that a certain class of AGNs has 
a higher Eddington ratio (ratio of bolometric to Eddington lu- 
minosity) among the AGN population; Narrow Line Seyfert 1 
galaxies (NLSls) and their high luminosity analogue, Narrow- 
Line QSOs (NLQs). They exhibit relatively narrow Balmer 
lines only slightly broader than the forbidden lines (FWHM of 
Hy8 < 2000 km s~'), optical Fe II multiplet emission, and low 
[OIII]-Hfi flux ratio, < 3 (see Pogge 2000). Strong soft X-ray 
excess (e.g. Pounds, Done & Osborne 1995; Boiler, Brandt & 
Fink 1996; Boiler et al. 2003), and rapid X-ray variability (e.g. 
Otani et al. 1996; Leighly 1999; Boiler et al. 2002; Gallo et al. 
2004) are often observed. Those characteristics indicate that 
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they have relatively small M Bll and larger M/(L E dd/c 2 ) (e.g. 
Brandt & Boiler 1998; Hayashida 2000; Mineshige et al. 2000). 

At super-Eddington accretion regime, the radiative effi- 
ciency per unit mass accretion is expected to decrease due to 
the onset of photon trapping (Begelman 1978). As a result, 
the emergent luminosity from an accretion flow starts to sat- 
urate at around a few times Lead (Abramowicz et al. 1988). 
Thus, rapid BH growth is possible in super-Eddington sources, 
NLSls/NLQs, (Rees 1992; Mineshige et al. 2000; Collin et al. 
2002; Kawaguchi 2003). However, the amount of BH growth 
via super-Eddington accretion is quite ambiguous, since it 
strongly depends on the assumed duration of NLSls/NLQs. 
The terms, "super-Eddington" and "sub-", are used in this 
Letter when an accretion rate of an AGN is larger and smaller 
than ~ 10LEdd/c 2 , where the emergent bolometric luminosity 
is about (0.5 - 1) LEdd- 

In this Letter, we assume that the ratio of the durations of a 
NLS1/NLQ and a BLSl/quasar may be proportional to the rel- 
ative number ratio of these two classes. Then, a simple algebra 
taking into account the fraction and average accretion rate of 
NLSls/NLQs is discussed, in order to evaluate qualitatively to 
what extent BHs can grow in super-Eddington accreting AGNs 
(§2). It is followed by several discussions in the final section. 



2. Super-Eddington accretion and BH growth 

2.1. Fraction and duration of NLS1 phase 

A fraction of AGNs are NLSls and NLQs: ~ 11% objects 
among a heterogeneous sample (Marziani et al. 2003), ~ 15% 
among the Early Data Release of the Sloan Digital Sky Survey 
(SDSS; Williams et al. 2002), and ~ (31 - 46)% among soft 
X-ray selected AGNs (Grupe et al. 1999; Grape 2004; Salvato, 
Greiner & Kuhlbrodt 2004). The fraction can be even higher 
in Extreme UV band; 8 NLSls out of 14 radio-quiet AGNs 
(Edelson et al. 1999). This wavelength dependency is likely 
due to the difference of the spectral energy distribution between 
NLSls/NLQs and BLSls/quasars (Boiler et al. 1996; Wang 
et al. 1996; Laor et al. 1997; see Fig. 14 of Kawaguchi 2003). 
Most of those NLSls have a redshift z smaller than 0.5. Here, 
we employ two numbers, 10% and 30%, as the relative fraction 
of the NLS1/NLQ among AGNs. 

In general, the duration (or the sum of multiple episodic 
phases) of an AGN is thought to be ~ 10 8 yr from several in- 
dependent arguments (e.g. Martini 2004; Yu & Tremaine 2002; 
Marconi et al. 2004). The duration of each episodic phase must 
be longer than 10 4 yr, to explain the proximity effect (e.g. 
Bajtlik, Duncan & Ostriker 1988) and the sizes of ionization- 
bounded narrow-line regions (Bennert et al. 2002). 

If statistical properties of AGNs (e.g. relative number ratio 
of super-Eddington to sub-Eddington sources) are quasi-steady, 
the fraction of NLS 1 s inferred above will be proportional to the 
duration of NLS1/NLQ phase. In fact, the evolution of AGNs 
is not so rapid in this redshift range, z = 0.5 — » 0, taking 5 Gyr: 
the number density decreases by a factor of 3, and the knee of 
LFs stays at almost the same luminosity (e.g. Boyle et al. 2000; 
Miyaji et al. 2001; Ueda et al. 2003). 



Then, let us simply assume that each NLS 1 (or NLQ) has a 
duration of ~ 10 7 yr or ~ 3 x 10 7 yr. This number (~ 10 7 yr) was 
briefly discussed together with the fraction of NLSls/NLQs 
also in a thesis (Grupe 1996). A characteristic timescale for 
BH growth can be roughly evaluated by (Salpeter 1964) 

M BH /M = 4.5 x 10 8 [j^2 ) y ear - (!) 

If BH growth mainly occures in NLS1/NLQ phase, the 
timescale of the order of 10 7 yr implies that M > 45 Z^Edd/c 2 . 

2.2. Accretion rate of NLSls 

The BH mass and gas accretion rate of AGNs onto their cen- 
tral BHs can be estimated using the observed relation between 
optical luminosity and the size of broad line region (Wandel, 
Peterson & Malkan 1999; Kaspi et al. 2000). Comparing opti- 
cal luminosity of NLSls and NLQs between observation and 
theoretical models, along with BH masses inferred from Hp 
and [O III] widths, their accretion rate is about 100L E dd/c 2 , on 
average (Fig. 1 1 in Kawaguchi 2003; Collin et al. 2002). 

Alternative and independent method to evaluate M Bii of 
NLSls/NLQs is spectral fitting of slim disk models at soft X- 
ray range (Wang & Netzer 2003). Moreover, if one tries to fit 
the broad-band spectra at optical-X-ray band, not only M B h 
but also M are estimated simultaneously (Kawaguchi 2003; 
Kawaguchi, Pierens & Hure 2004; Kawaguchi, Matsumoto & 
Leighly in prep.; see Kawaguchi 2004). Currently, these mod- 
els are the only ones that explain broadband SEDs of NLSls, 
in our best knowledge. This kind of estimations indicate Mbh 
and M which are similar to those derived by the line widths in 
the previous paragraph. Spectral fitting is relatively time con- 
suming, but this success implies that we can use the line widths 
for quick estimations, as the first order. 

2.3. Linear growth of BHs 

Now, we start with an assumption that the gas accretion rate 
onto a central BH is almost constant: M « const. For instance, 
a following scenario is possible for the constant-M case: large 
amount of gas accretion (super-Eddington accretion phase = 
NLSl/NLQ-phase) may start induced by an unknown trigger. 
With a growing BH mass, the accretion rate in unit of the 
Eddington rate will be reduced down to sub-Eddington phase 
(broad line Seyferts and quasars). Accretion rate (in Eddington 
unit) will decrease continuously, and eventually it will become 
lower than ~ 0.2L E dd/c 2 , followed by a transition from a radia- 
tively efficient to an inefficient flow (i.e. optically-thin advec- 
tive flow; e.g. Meyer, Liu & Meyer-Hofmeister 2000), such as 
Low Luminosity AGNs. 

Employing this view, the growth of Mbh with time t is lin- 
ear as is expressed below: 

M BH « = M BH (t = 0) + M X t. (2) 

Substituting the average accretion rate of NLSls/NLQs to M, 
we have 

M BH (t = 10 7 yr) * 8 x M BH (t = 0), and (3) 
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M BH (f = 3 x 10 7 yr) * 23 x M BH (t = 0), (4) 
ifM^100 LEdd( ^ 2 106 - 5yr) . 

On the other hand, in a sub-Eddington phase the factor of 
BH growth during ~ 10 8 yr is not so significant: 



M BH (f = 10 8 yr) = 1.8 x M BH (r = 0), 



if M = 3 



= 10 75 yr) 



2.4. Exponential growth of BHs 

Next, another assumption on accretion rate is that M /(^Edd/c 2 ) 
may be constant with time by some unknown self-regulation 
mechanism(s). In this case, sub-Eddington phase does not fol- 
low at the end of super-Eddington phase, contrary to the previ- 
ous case. Instead, in each tagger of AGN activity, M /(^Edd/c 2 ) 
is (somehow randomly) chosen: if a nuclei has a super- 
Eddington (sub-) accretion rate it will appear as a NLS1/NLQ 
(a broad-line Seyfert/quasar). The growth of M B h will be ex- 
ponential with this assumption, and is written as, 



M BH (0 = M BH (f = 0) x exp 
= M BH (f = 0) x exp 



t x 



M 
M BH 
t 



(. 



M 



4.5 x 10 8 yr/UEdd/c 2 
Substituting the average accretion rate (§2.2), we get 
M BH (r = 10 7 yr) * 9 x M BR (t = 0), and 
M BH (f = 3 x 10 7 yr) * 800 x M BH (f = 0). 



(5) 

(6) 
(7) 



Given the average accretion rate (~ 3LEdd/c 2 ) and duration 
(~ 10 8 yr) of BLSls/quasars, a BH can grow by a factor of 1.9 
in the sub-Eddington phase. 

3. Discussions and conclusions 

In principle, a BH mass based on line widths is the mass in- 
side the broad line region, and does not necessarily equal to the 
true Mbh ■ Outer part of accretion disks in NLS 1 s/NLQs, where 
optical continuum is emitted, are self-gravitating (Kawaguchi 
et al. 2004). Such disks are quite massive, and M B h based on 
line widths can be systematically overestimated by a factor of 
~ 2. If this is taken account into M B h and M estimations, the 
average accretion rate of NLS 1 s/NLQs might be 200LEdd/c 2 , 
rather than 100LEdd/c 2 . Thus, the factor of BH growth pre- 
sented in §2.4 (eqs. [6J will be multiplied by themselves. 
While, it is too efficient to grow a BH with this accretion rate 
for a constant-M case (§2.3). For sub-Eddington systems, on 
the other hand, there is no constraints upon the presence of 
outer, massive, self-gravitating disks, since such region (if ex- 
ists) radiates at Near-IR where the emission from dusty torus 
dominates. 

We have utilized the relative fraction and the average accre- 
tion rate of NLSls/NLQs obtained for z < 0.5. Apparently, the 
Eddington ratios of AGNs detected in SDSS tend to increase 



with redshift (McLure & Dunlop 2003). Thus, BH growth via 
super-Eddington accretion would be much more important in 
the young universe than discussed above. 

The tight correlation between spheroid mass of galaxies 
and BH mass indicates that they co-evolve. Indeed, high star 
formation activities in host galaxies of NLSls/NLQs are in- 
ferred (e.g. Moran, Halpern & Helfand 1996; Kawaguchi & 
Aoki 2001). Since the gas accretion rate of NLSls/NLQs is 
~ 2M yr~' at most (Collin & Kawaguchi 2004), a small frac- 
tion of mass loss by starburst streaming toward the central BH 
is enough to form a NLS1/NLQ. 

We note that the super-Eddington accretion phase (with 
L < ^Edd) discussed here is not a surprising assumption in 
the context of cosmic structure formation history. In semi- 
analytical approaches (e.g. Kauffmann & Haehnelt 2000), for 
instance, most of accretion events happen in super-Eddington 
regime if the episodic phase of an AGN is assumed to be less 
than 10 8 yr (Cattaneo 2001). 

For quasars, with their bolometric luminosity L > 10 46 
erg/s, LFs can be obtained by direct observations to high red- 
shifts in a wide range of wavebands. Comparing the integra- 
tion of those LFs with the volume mass density of BHs in the 
local universe, the radiative efficiency seems to be quite high 
(> 0.15) in those luminous objects (e.g. Yu & Tremaine 2002; 
Elvis, Risaliti & Zamorani 2002). Such a high efficiency indi- 
cates that most BHs of quasars are rotating with sub-Eddington 
accretion rates. However, the mass density of local BHs are 
dominated by quite massive BHs at the knee of mass functions 
or LFs. Therefore, constraints upon BH growth based on above 
discussions are likely valid at high mass end; e.g. BH growth 
from 10 8 M Q to 1O 1O M . 

On the other hand, BH growth for lower luminous objects 
with M B h < 1O 8 M is quite unclear. The way BHs grow in 
Seyferts is not necesarilly the same as that of quasars. Actually, 
Eddington ratios of AGNs seems to be higher for objects with 
smaller M B h (e.g. Collin et al. 2002; Collin & Kawaguchi 
2004). Yu & Tremaine (2002) argued that relatively low lumi- 
nous AGNs should have less radiative efficiency, < 0.1. Thus, 
super-Eddington accretion (which is expected to have very low 
efficiency) seems to play an important role in BH growth from 
seed BHs to ~ 1O 8 M . This is not in conflict with the idea 
above that BHs in quasars grow with high radiative efficiency 
(via sub-Eddington accretion). 

Given that a BH mass increases by the factors mentioned 
above, is it enough to grow BHs found in quasars with ~ 
1O 8 ~ 1O M from seed BHs? It may be possible to make such 
supermassive BHs if i) the mass of seed BHs at the center of 
galaxies is large enough (e.g. > 10 M ), and/or ii) the duration 
and average accretion rate of super-Eddinton accreting AGNs 
are much larger at high redshift than those estimated for local 
objects. Instead, iii) BH growth at the highest mass end may 
be indeed dominated by merging events as expected in semi- 
analytical approaches (Cattaneo 2002). 

In summary, we have argued that the growth of black hole 
(BH) mass by factor of 8 - 800 happens in the super-Eddington 
accretion phase. There is a possible, systematic underestima- 
tion of M /(Lsdd/c 2 ) by a factor of ~ 2 due to an overes- 
timation of BH mass by massive accretion discs for super- 
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Eddington objects. If it is true, the factor of BH growth above 
may be larger by order(s) of magnitude. On the other hand, 
the growth factor expected in sub-Eddington active galactic nu- 
clei (AGNs) is only ~ 2. Therefore, the cosmic BH growth 
by accretion is dominated by super-Eddington phase, rather 
than sub-Eddington phase which is the majority among AGNs. 
These values are based on the fraction and the average accre- 
tion rate of Narrow-Line Seyfert 1 galaxies and Narrow-Line 
quasars obtained for z < 0.5. If those numbers are larger at 
higher redshift (where BHs are less grown), BH growth via 
super-Eddington accretion would be even more important. 

For further analysis, we need i) volume limited samples of 
AGNs to derive an accurate relative fraction of NLSls/NLQs, 
and ii) constraints on z-dependency of the relative fraction and 
the average M/(LEdd/c 2 ) to determine the importance of super- 
Eddington accretion in BH growth history. In this Letter, we 
divide radio-quiet type I AGNs into two subclasses: objects 
with M > or < lOL E dd/c 2 - It would be more appropriate to 
divide into fine bins for dealing with the observed continuous 
distribution of M/(L E ^/c 2 ); e.g. bins for 10 - 100L E dd/c 2 and 
100 - lOOOLEdd/c 2 ' etc. For this a i m > Si) more sophysticated 
(with less uncertainty) methods for M/(L E ^/c 2 ) estimations 
are necesary. 
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